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Abstract 

Exponentiating the hypergeometric series qFl(1, 1, . . . , 1; z), L = 0, 1, 2, ... , 
furnishes a recursion relation for the members of integer sequences &z,(n), n = 
0, 1, 2, . . .. For L > 0, the &k(n)'s are certain generalizations of conventional 
Bell numbers, bo(n). The corresponding associated Stirling numbers of the 
second kind are also generated and investigated. For L = 1 one can give 
a combinatorial interpretation of the numbers b\{n), and of some Stirling 
numbers associated with them. We also consider the L > 1 analogues of Bell 
numbers for restricted partitions. 

The conventional Bell numbers [|l[] b (n), n = 0,1, 2, .. ., have a well known exponential 
generating function 

oo n 

5oW = e( eZ - 1 )=^6o(n)^ [ , (1) 

n=0 

which can be derived by interpreting bo(n) as the number of partitions of a set of n dis- 
tinct elements. In this note we obtain recursion relations for related sequences of positive 
integers, called 6i(n), L = 0, 1, 2, . . . , obtained by exponentiating the hypergeometric series 
oF L (l,l,...,l^) defined by @: 
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(2) 



(for which we shall use throughout the short notation Fl(z)) and which includes the special 
cases o-^q(^) = e * an d oFi(z) = ^o(2\/S), where Io(x) is the modified Bessel function of 
the first kind. For L > 1, the functions oFl(z) are related to the so-called hyper-Bessel 
functions ||, [|], 0, which have recently found application in quantum mechanics @, [0- 
Thus, we are interested in b^n) given by 



e [ F L (z)-l] 



EWt^TT' (3) 



n 

n=0 



thereby defining a hypergeometric generating function for the numbers &L(n). From eq.([3D 
it follows formally that 

dz r ' 

For L = the r.h.s of eq.(HI) can be evaluated in closed form: 



Mn) = (n!) £ .^( e ^W- 1 l)U. (4) 



k=0 



z— I e 



(5) 



The first equality in is the celebrated Dobihski formula [p]], ||, f9|. The second equality 
in eq.(^[) follows from observing that for a power series R(z) = Yl'kLo ^k zk 

A \ n 00 



k=0 



holds, and by applying eq.(^) to the exponential series (Ak = (A;!) -1 ). 

The reason for including the divisors (n\) L+1 rather than n\ as in the usual exponential 
generating function, arises from the fact that only through eq.(|J) are the numbers 6t(n) 
actually integers. This can be seen from general formulas for exponentiation of a power 
series |J, which employ the (exponential) Bell polynomials, complicated and rather unwieldy 
objects. It cannot however be considered as a proof that the biiji) are integers. At this 
stage we shall use eq.(|3]) with biin) real and apply to it an efficient method, exposed in 
||, which will yield the recursion relation for the b^n). (For the proof that the b^n) are 
integers, see below eq. flTT|) ) ■ To this end we first obtain a result for the multiplication of two 
power-series of the type (^|). Suppose that we have to multiply f(x) = J2™=o a L( n ) ^l+i 
and g(x) = J]~ =0 c L (n) py£+r- We get f{x) ■ g{x) = J2n=o d Mj^fi+T, where 



= £ M^ = ±{;) L+1 a4r)c L (n-r). (7) 

r+s=n ^ ' ^ ' r=0 
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Substitute eq.([|) into eq.(|3]) and take the logarithm of both sides of eq.( 



n=l 



v n=0 



(n!) 



£+1 



(8) 



Now differentiate both sides of eq.(H) an d multiply by z. It produces 



v n=0 



i£+l 



E 

v n=0 



n 



n!) 



L+l 



n=0 



which with eq.(^) yields the desired recurrence relation 



b L (n + l) 



n 

— Y 



n _l_ ]_ ^ ( ft 

fc=0 



L+l 



E 



fc=0 



n \ / n + 1 



(n+l-k) b L {k), 



b L (k), 



n!) 



n = 0,1,... 



(9) 



6l(0) = 1. 



(10) 

(11) 
(12) 



Since eq. ([[]]) involves only positive integers, it follows that the &L(n) are indeed positive 
integers. For L = one gets the known recurrence relation for the Bell numbers H: 



n 

b (n + l) = Y,(l)bo(k). 



(13) 



fc=0 



We have used eq.flnp to calculate some of the 6^(n)'s, listed in Table I, for L — 0, 1, . . . , 6. 
Eq. (|TTD , for n fixed, gives closed form expressions for the fox(n) directly as a function of L 
(columns in Table I): b L {2) = 1 + 2 L , 6 L (3) = 1 + 3 • 3 L + (3!) L , 6 L (4) = 1 + 4 • 4 L + 3 • 6 L + 
6 • 12 L + (4!) L , etc. 

The sets of b^n) have been checked against the most complete source of integer sequences 
available [ 10| . Apart from the case L = (conventional Bell numbers) only the first non- 
trivial sequence L — 1 is listed: it turns out that this sequence b\ (n), listed under the 
heading A023998 in [|10], can be given a combinatorial interpretation as the number of block 
permutations on a set of n objects, which are uniform, i.e. corresponding blocks have the 
same size |fT2]|. 

Eq. (jl|) can be generalized by including an additional variable x, which will result in 
"smearing out" the conventional Bell numbers bo(n) with a set of integers So(n,k), such 
that for k > n, S (n, k) = 0, and S (0, 0) = 1, S (n, 0) = 0. In particular, 



B (z,x) = e 



x(e 



L) = E 



n=0 



x 



k=l 



n! 



(14) 



which leads to the (exponential) generating function of So(n,l), the conventional Stirling 
numbers of the second kind, (see IfEfl, HI), in the form 



3 



1)' 



E 

n=l 



S (n,l) 



nl 



(15) 



and defines the so-called exponential or Touchard polynomials In (x) as 

n 

l n 0) (x) = Y,So(n,k)x k . 

k=l 

They satisfy 

4 0) (i) = ^oH, 



(16) 



(17) 



justifying the term "smearing out" used above. 

The appearance of integers in eq.(|3|) suggests a natural extension with an additional 
variable x: 



B L (z,x) = e x[oFL(z) - 1] = J2 



n=0 



^2S L (n, k) 



X 



k=l 



(n!) 



L+l 



(18) 



where we include the right divisors (n!) L+1 in the r.h.s of fljjf). 

This in turn defines "hypergeo metric" polynomials of type L and order n through 



x 



(19) 



k=l 



which satisfy 



J?>(1) = &l(ti), 



(20) 



with the &l(™) of eg- fllCf ). Thus, the polynomials of eq. ([l9|) "smear out" the 6l("-) with 
the generalized Stirling numbers of the second kind, of type L, denoted by Sl(ti, k) (with 
Sl(ji, k) — 0, if k > n, Sl{ti, 0) = if n > and Sl(0, 0) = 1), which have, from eq . (|18D the 
"hypergeo metric" generating function 



MM - iy 



E 

n=l 



Sl(u,1) 



L+l 



L = 0,l,2,.... 



(21) 



Eq.pT|) can be used to derive a recursion relation for the numbers Sl(u, k), in the same 
manner as eq.([|) yielded eq.(|12|). Thus we take the logarithm of both sides of eq. (pl|), 
differentiate with respect to z, multiply by z and obtain: 



E 

v n=0 



S L {n,l- 1) 



\L+1 



E 

v n=0 



7V. 



E 

n=0 



nS L (n,l) n 
(n\) L+1 Z ' 



(22) 
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which, with the help of eq.(|7]), produces the required recursion relation 
S L (n + l,l)= (I) (T) S L (k,l-l), 



k=l-l 

^(0,0) = !, S L (n,0) = 0, 



(23) 
(24) 



which for L = is the recursion relation for the conventional Stirling numbers of the second 
kind UU, ||, and in eq.( p3|) the appropriate summation range has been inserted. Since the 
recursions of eq.(^) and eq. (|24"l) involve only integers we conclude that Sl(ti, I) are positive 
integers. 

We have calculated some of the numbers Si(n,l) using eq . fl2l"D and have listed them in 
Tables II and III, for L = 1 and L = 2 respectively. Observe that Si(n,2) = ^ ^^-i ) — ^ 
and SL(n,n) = {n\) L , L = 1,2. Also, by fixing n and Z, the individual values of Sl(ti,1) 
have been calculated as a function of L with the help of eq. (p3[) , see Table IV, from which 
we observe 



S L (n,n) = (n\y 



1,2, 



(25) 



which is the lowest diagonal in Table IV. We now demonstrate that the repetitive use of 



eq.(|23]) permits one to establish closed- form expressions for any supra-diagonal of order p, i.e. 
the sequence Sl(ti +p,n), for p — 1, 2, 3, . . ., if one knows the expression for all Sl(ti + k, n) 
with k < p. We shall illustrate it here for p — 1,2. To this end fix I = n on both sides 
of eq.(p3|). It becomes, upon using eq . (|25|) , and defining Ol(ti) = Sl(ti + l,n), a linear 
recursion relation 



a L ln 



with the solution 



n[(n + 1)! 



+ (n+ l) L a L (n- 1), 



«l(0) = 0, 



a L (n) = S L (n + l,n) 



n(n + 1) 
2 

(n + l)\ 



(n + 1) 



S (n + l,n), 



(26) 

(27) 
(28) 



which gives the second lowest diagonal in Table IV. Observe that for any L, Sl(ji + 
l,n) is proportional to S (n + l,n) = n{n + l)/2. The sequence S\{n + l,n) = 
1, 9, 72, 600, 5 400, 8 564 480, ... is of particular interest: it represents the sum of in- 
version numbers of all permutations on n letters ||10|| . For more information about this and 



related sequences see the entry A001809 in flQfl . The Sl(u + 1, n) for L > 1 do not appear to 
have a simple combinatorial interpretation. A recurrence equation for /3i(n) = Sl(to + 2, n) 
is obtained upon substituting eq.(^Bl) and eq.(^7|) into eq.(^3|): 
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n(n + 1) 
2! 



(w + 2) 
2! 



to- 1 1 



2 L 3 L 



-r +(to + 2)^(to-1), /3 l (0) = 0. (29) 



It has the solution 



S L (n + 2,n) 



n(n + l)(n + 2) 
3^23 



(to + 2) 



3 / ,N 4 

-r(TO - 1) + — 

2 L v ; 3 L 



(30) 



which is a closed form expression for the second lowest diagonal in Table IV. Clearly, eq.(pO| 
for L = gives the combinatorial form for the series of conventional Stirling numbers 



5*0 (to + 2, to) 
In a similar way we obtain 



to(to+ 1)(to + 2)(3to + 1) 



S L (n + 3, to) 



to(to + 1)(to + 2) (to + 3) 



x n 



3 • 2 4 

L 



+ TO 



4! 



(n + 3)! 
3 

3 L+1 



8 L 



+ 



2 + 2 • 3 i 



8 L 



which for L = reduces to 



(31) 



(32) 



So(to + 3,to) = — to 2 (to + 1) 2 (to + 2)(to + 3). 
48 



Combined with the standard definition 



So W) = ^ £(-D l (I)*' 

fc=l 



(33) 



(34) 



eqs.( p8|) , ( ^T|) and (|33D give compact expressions for the summation form of So (to + P, to). 
Further, from eq.(|31), use of eq.(H) gives the following generating formula 



S Q (n,l) 



V 



/! 



d 

dz 

4) [(l-=)'-H 



,fc=l 



2=1 



2=1 
TO > /. 



(35) 
(36) 



The formula (|I|) can be generalized by putting restrictions on the type of resulting par- 
titions. The generating function for the number of partitions of a set of to distinct elements 



without singleton blocks bo(l,n) is ||, |T4|, [15 



~7i 

B (l,z) = e eZ - 1 - z = J2bo(^n)-, 



n=0 



(37) 



or more generally, without singleton, doubleton . . ., p— blocks (p — 0, 1, . . .) is [[yj 



6 



J. 



B (p, z) = e eZ " s "=o *r = ^ b (p, n) — , (38) 



n=0 



with the corresponding associated Stirling numbers defined by analogy with eq.(|l4|) and 
eq. fl22|) . The numbers b (l,n), b (2,n), b (3,n), 6 (4, n) can be read off from the sequences 
A000296 , A006505 , A057837 and A057814 in 0, respectively. For more properties of these 



numbers see 11 



We carry over this type of extension to eq.(|) and define bi{p,n) through 

k < ~*~ > n 

B L (p,z) = e °^«-n=o<sifcr =J2h(p,n)— * (39) 

n=0 \ '> 

where &r,(0, n) = bi(n) from eq.(f|). (We know of no combinatorial meaning of &z\(p,n) for 
L > 1, p > 0). The &i(p, n) satisfy the following recursion relations: 



n-p L 

fc(p>") = £(;)(T) b L (p,h), (40) 



fc=0 



& L (p,0) = l, (41) 
b L (p, 1) = 6 x (p, 2) = ■ • ■ = b L (p,p) = 0, (42) 
6 £ (p,p + l) = l. (43) 

That the &c(p, n) are integers follows from eq.fl4"0|). Through eq.(|3"9"|) additional families of 
integer Stirling-like numbers Si tP (n, k) can be readily defined and investigated. 

The numbers bo(p, n) are collected in Table V, and Tables VI and VII contain the lowest 
values of bi(p,n) and 62 (Pj n), respectively. 

Formula ([3]) can be used to express e in terms of bo(n) in various ways. Two such lowest 
order (in differentiation) forms are 

e=1+ J±m) = (44) 



m 

v n=0 



00 



= i n(E M^ti) . (45) 

\n=0 ' ' / 

In the very same way, eq.(^) can be used to express the values of qFl{z) and its derivatives 
at z = 1 in terms of certain series of 6^(n)'s. For L = 1, the analogues of eq.(f44"D and eq. (f4~5l) 
are 

7 (2) = l + lng^, (46) 

OO 



Jo(2) + ln(J 1 (2)) = l + ln(g (K 6l |" 1 ^ ! ) )2 ) (47) 
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and for L = 2 the corresponding formulas are 



) F 2 (l,l;l) = l + ln [J2 



b 2 (n) 



v n=0 

oo 



n! 



)F 2 (1, 1; 1) + In ( F 2 (2, 2; 1)) = 1 + In ^ 



v n=0 



& 2 (^ + 1) 

n + l) 2 (n!) 3 



(48) 
(49) 



By fixing 2:0 at values other than Zq = 1, one can link the numerical values of certain 
combinations of o-^i(l, 1, • • • ; ^o) > o-^l(2, 2, . . .; zq),. . . and their logarithms, with other series 
containing the 6i(ri)'s. 

The above considerations can be extended to the exponentiation of the more general 
hypergeometric functions of type F L (ki, k 2 , . . . , k^, z) where kx, k 2 , ■ ■ ■ , &l are positive in- 
tegers. We conjecture that for every set of fc n 's a different set of integers will be generated 
through an appropriate adaptation of eq.(§). We quote one simple example of such a series. 
For 



oF 2 (l,2;z) = ]T 



n=0 



[n + l)(n\) 3 



(50) 



eq.(||) extends to 



e [ F 2 (l,2;z)-l] = ^f 2i 



n 



n=0 



(51) 



where 



f 2 (n) = (n + l)(n\f 



d n 
dz n 



J F 2 (1,2; Z )-1] 



(52) 



2 = 



turn out to be integers: / 2 (n), n = 0, 1, . . . , 8 are: 1, 1, 4, 37, 641, 18 276, 789 377, 48 681 
Oil, etc. The analogue of equations (|23|) and ( flip is: 



F 2 (l,2;l) = l + lnPr 



v n=0 



[n + l)(n!) 3 



(53) 
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TABLES 



TABLE I. Table of b L (n); L, n = 0, 1, . . . , 6. 



L 


Mo) 


Mi) 


M2) 


M3) 


M 4 ) 


M5) 


M6) 





i 


i 


2 


5 


15 


52 


203 


1 


i 


i 


3 


16 


131 


1 496 


22 482 


2 


i 


i 


5 


64 


1 613 


69 026 


4 566 992 


3 


i 


i 


9 


298 


25 097 


4 383 626 


1 394 519 922 


4 


i 


i 


17 


1 540 


461 105 


350 813 126 


573 843 627 152 


5 


i 


i 


33 


8 506 


9 483 041 


33 056 715 626 


293 327 384 637 282 


6 


i 


i 


65 


48 844 


209 175 233 


3 464 129 078 126 


173 566 857 025 139 31 







TABLE II. Table 


of S L (n,l); 


L = 1, l,n 


= 1,2,... 


,8. 




I 


5i(l,0 


Si (2,0 


Si (3,0 


Si (4,0 


Si (5,0 


Si (6,0 


Si (7,0 


Si (8,0 


1 


1 


1 


1 


1 


1 


1 


1 


1 


2 




2 


9 


34 


125 


461 


1 715 


6 434 


3 






6 


72 


650 


5 400 


43 757 


353 192 


4 








24 


600 


10 500 


161 700 


2 361 016 


5 










120 


5 400 


161 700 


4 116 000 


6 












720 


52 920 


2 493 120 


7 

8 














5 040 


564 480 
40 320 



TABLE III. Table of S L (n, l); L = 2, I, n = 1, 2, . . . , 8. 
I 5 2 (1,0 S 2 (2,0 S 2 (3,0 S 2 (4,0 S 2 (5,0 S 2 (6,0 S 2 (7,0 S 2 (8,0 



1111111 1 1 

2 4 27 172 1 125 7 591 52 479 369 580 

3 36 864 17 500 351 000 7 197 169 151 633 440 

4 576 36 000 1 746 000 80 262 000 3 691 514 176 

5 14 400 1 944 000 191 394 000 17 188 416 000 

6 518 400 133 358 400 23 866 214 400 

7 25 401 600 11 379 916 800 

8 1 625 702 400 
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TABLE IV. Table of S L (n,l); I, n = 1, 2, . . . , 6. 



I S L (l,l) S L (2,l) S L (3,l) S L (4,l) S L (5,l) S L (6,l) 



11111 1 1 

2 (2!) L 3 • 3 L 4 • 4 L + 3 • 6 L 5 • 5 L + 10 • 10 L 6 • 6 L + 15 • 15 L + 10 • 20 L 

3 (3!) L 6 • 12 L 10 • 20 L +15 • 30 L 15 • 30 L + 60 • 60 L + 15 • 90 L 

4 (4!) L 10 • 60 L 20 • 120 L + 45 • 180 L 

5 (5!) L 15 • 360 L 

6 (6!) L 



TABLE V. Table of b (p, n); p = 0, 1, 2, 3; n = 0, . . . , 10. 



n 


6o(0,n) 


6 (l,n) 


&o(2,n) 


&o(3,n) 





1 


1 


1 


1 


1 


1 











2 


2 


1 








3 


5 


1 


1 





4 


15 


4 


1 


1 


5 


52 


11 


1 


1 


6 


203 


41 


11 


1 


7 


877 


162 


36 


1 


8 


4 140 


715 


92 


36 


9 


21 147 


3 425 


491 


127 


10 


115 975 


17 722 


2 557 


337 
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TABLE VI. Table of b^p, n); p = 


0,1,2; n = 0,...,9. 




n 


6i(0,n) 


&i(l,n) 


bi (2, n) 





1 


1 


1 


1 


1 








2 


3 


1 





3 


16 


1 


1 


4 


131 


19 


1 


5 


1 496 


101 


1 


6 


22 482 


1 776 


201 


7 


426 833 


23 717 


1 226 


e 

o 




^1 *i Q71 


00 1 


9 


277 006 192 


11 893 597 


493 333 




TABLE VII. Table of b 2 (p, n); p = 


= 0,1,2; n = 0,...,8. 




n 


b 2 {0,n) 


62(1, n) 


62(2, n) 





1 


1 


1 


1 


1 








2 


5 


1 





3 


64 


1 


1 


4 


1 613 


109 


1 


5 


69 026 


1 001 


1 


6 


4 566 992 


128 876 


4 001 


7 


437 665 649 


4 682 637 


42 876 


8 


57 903 766 800 


792 013 069 


347 117 
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